The information content of a new observable: the case of the nuclear neutron skin 
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We address two questions pertaining to the uniqueness and usefulness of a new observable: (i) 
Considering the current theoretical knowledge, what novel information does new measurement bring 
in? (ii) How can new data reduce uncertainties of current theoretical models? We illustrate these 
points by studying the radius of the neutron distribution of a heavy nucleus, a quantity related to 
the equation of state for neutron matter that determines properties of nuclei and neutron stars. By 
systematically varying parameters of two theoretical models and studying the resulting confidence 
ellipsoid, we quantify the relationships between the neutron skin and various properties of finite 
nuclei and infinite nuclear matter. Using the covariance analysis, we identify observables and pseudo- 
observables that correlate, and do not correlate, with the neutron skin. By adding the information on 
the neutron radius to the pool of observables determining the energy functional, we show how precise 
experimental determination of the neutron radius in 208 Pb would reduce theoretical uncertainties 
on the neutron matter equation of state. 

PACS numbers: 21.10.Gv, 21.60.Jz, 21.65.Mn 25.30.Bf 
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Introduction. Nuclei communicate with us through a 
great variety of observables. Some are easy to measure, 
some take a considerable effort and experimental inge- 
nuity. Without any preconceived knowledge, all differ- 
ent observables are independent of each other and can 
usefully inform theory. On the other extreme, new data 
would be redundant if our theoretical model were perfect. 
Reality lies in between. In this paper we show how to as- 
sess the uniqueness and usefulness of an observable, i.e., 
its information content with respect to current theoreti- 
cal models. We also quantify the meaning of a correlation 
between different observables. 

Let us consider a model characterized by a number 
of parameters p = (pi, ...,pf) defining the model space. 
Those parameters may be, e.g., coupling constants of the 
effective Hamiltonian and effective charges characteriz- 
ing operators in the assumed Hilbcrt space. Calculated 
observables are functions of these parameters. Since the 
number of parameters is usually much smaller than the 
number of observables, correlations exist between com- 
puted quantities. Moreover, since the model space has 
been optimized to a limited set of observables, there may 
also exist correlations between parameters. 

The paper is organized as follows. We first explain 
the statistical approach used to estimate theoretical un- 
certainties and quantify correlations between calculated 
observables. We next discuss the importance of the neu- 
tron skin measurement, briefly review related theoretical 
work, and describe the model used and the set of fit- 
observables employed. The remainder of the paper con- 
tains results, conclusions, and the outlook for the future. 

Regression analysis. Usually, most of the model space 
produces observables which are far from reality. There- 
fore, one needs to confine the model space to a "phys- 



ically reasonable" domain. That can be achieved by a 
least squares regression analysis. To this end, one selects 
a pool of fit-observables O which are used to calibrate 
p. The optimum parameterization po is determined by 
a least squares fit with the global quality measure 
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where "th" stands for the calculated values, "exp" for ex- 
perimental data, and AO for adopted errors. Having de- 
termined po , an expectation value of an observable A can 
be computed at A(po). However, there remain uncertain- 
ties, originating both from the errors in fit-observables 
and from a limited reliability of the model. To estimate 
the root-mean-square (rms) variation of A, one needs to 
define a physically reasonable domain around po- Near 
the minimum, the x 2 -landscape is given by a confidence 
ellipsoid (see Sec. 9.8 of Q): 
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The physically reasonable domain p is defined as that 
multitude of parameters around po that fall inside the 
covariance ellipsoid x 2 — \n + 1, i.e., 



(p-po)Ai(p-Po) < 1. 
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In terms of Aiij, the covariance between two observables 
A and B becomes: 
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For A—B, Eq. ([5]) gives variance A 2 A that defines an 
uncertainty of an observable. In addition, one can also 
establish the Pearson product-moment correlation coef- 
ficient between two observables [ljj: 



cab = 



\AA AB\ 
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A value cab = 1 means fully correlated and cab = is 
totally uncorrelated. 

Nuclear neutron skin. The covariance analysis is a 
standard statistical tool that can be applied to any the- 
oretical model that has been optimized to the data. In 
this work, we illustrate the general concept by consider- 
ing one particular observable of fundamental importance 
for nuclear physics and astrophysics: the rms radius 
of the neutron density distribution in a heavy nucleus, 
r™ = (r 2 )l/ 2 . The size of r™ s is strongly correlated 
with many properties characterizing neutron-rich matter 
found in neutron-rich nuclei [2j and in neutron stars |3|. 
The highly anticipated Lead Radius Experiment (PREX) 
at Jefferson Laboratory will use the parity- violating elec- 
troweak asymmetry in the elastic scattering of polarized 
electrons to determine the neutron radius of 208 Pb with 
a projected experimental precision of 1%, in a model- 
independent fashion Q. Below, we apply the covariance 
analysis to address questions pertaining to neutron-rich 
matter in general and PREX experiment in particular: 
(i) What quantities that are experimentally accessible 
from finite nuclei correlate best, or do not correlate, with 
neutron radius? (ii) How robust are correlations between 
observables from finite nuclei and nuclear matter prop- 
erties (NMP)? (hi) To what extent would precise data 
on the neutron radius in 208 Pb enhance the predictive 
ability of theory? 

A quantity that is related to r™ s is the neutron skin 
r s kin = — r™ s @, Hj]. The usefulness of neutron 
skin lies in its strong dependence on the isovector den- 
sity pi — p n — p p and a much weaker dependence on 
the isoscalar, or total, density p = p n + p p . A num- 
ber of relationships, or correlations, have been estab- 
lished between r s ki n in heavy nuclei and various NMP 
and observables in finite nuclei (see @ for an early discus- 
sion). Those include: symmetry energy at the saturation 
point a sym (p e q) [E 0413) slope of bulk symmetry energy 
a sym = da sym /dp (proportional to the pressure differ- 
ence between neutrons and protons) at p cq 0, EH E3 
and at p=0.1 nucleons/fm 3 [hi]], slope of binding energy 
of neutron matter d(E/A) n /dp n at p n =0.1 neutrons/fm 3 
[Til [l5l - [T7j (proportional to the neutron pres- 
sure), the symmetry correction to the incompressibility 
AK Q, low-energy electric dipole strength attributed 
to the Pygmy Dipole Resonance (PDR) [18l - l20f . neu- 
tron form factor [8], and r^m in different nuclei [2lj . It 
has also been found that there are NMP which correlate 
poorly with r s ki n : equilibrium nuclear matter binding en- 
ergyand saturation density p eq @, incompressibility K 
[1, 11 1, and enhancement factor of the Thomas- Reiche- 



Kuhn sum rule ktrk (related to the isovector effective 
mass) 0- 

Some of the previous theoretical papers dealing with 
neutron skin correlations have explored the dependence 
between observables by explicit variation of selected 
properties (e.g. symmetry energy) within the given 
model, see, e.g., Refs. [1, E3- The present covariance 
analysis is the least biased and most exhausting way to 
find out the correlations ([5]) between all conceivable ob- 
servables. There remain, however, what is called system- 
atic errors which are here hidden constraints and limi- 
tations of the given model. Such systematic errors can 
only be determined by comparing different models or suf- 
ficiently flexible variants of a model. A comparison of 
different models as in, e.g., Refs. |, EE El is thus an 
instructive complement. But the use of different models 
is not appropriate to quantitatively assess the correlation 
between observables. For that reason, our study is based 
on the covariance analysis within the framework of one 
model. 

The Model. The theoretical approach employed in this 
study is the self-consistent mean field theory in the nu- 
clear Density Functional Theory (DFT) formulation [22j ■ 
At its heart lies the nuclear energy density functional 
(EDF) that is built from the nucleonic intrinsic densities 
and - in a relativistic version - meson fields. The nuclear 
DFT framework has been successful in describing a broad 
range of nuclear properties, including ground-state prop- 
erties, excited states, particle decays, and fission. Over 
the last few years, however, it has become evident that 
the standard functionals are too restrictive when one is 
aiming at the detailed quantitative description and ex- 
trapolability. Consequently, various strategies have been 
devised to develop realistic EDF of spectroscopic qual- 
ity [23| • Early attempts to employ statistical methods of 
linear-regression and error analysis [24| have been revived 
recently and been applied to determine the independence 
of EDF parameters, their errors, and the errors of calcu- 
lated observables. |13l |25I428| |. The major uncertainty in 
EDF lies in the isovector channels that are poorly con- 
strained by experiment. In this context, neutron skin 
data are crucial. 

The EDF used in this work is the Skyrme functional 
SV-min of Ref. It is characterized by F=14 cou- 

pling constants (listed in Table V therein). The observ- 
ables chosen to define % 2 during optimization of SV-min 
embrace nuclear bulk properties (binding energies, sur- 
face thicknesses, charge radii, spin-orbit splittings, and 
pairing gaps) for selected semi-magic nuclei which are 
proven to allow a reasonable DFT description. For a list 
of chosen observables, pseudo-observables, and adopted 
errors, see Tables I-IX of Ref. [H]. NMP are not in- 
cluded in the fit data for SV-min. This allows to count 
them as extrapolated observables in the present correla- 
tion study. The parameter set po of SV-min provides a 
very reasonable description of finite nuclei and nuclear 
matter (if=222MeV, a syBl =30.7MeV, effective nucleon 
mass m*/m=0.95). 
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A second EDF used in this study is the relativistic 
mean field. We use it here in a traditional form in 
which Dirac nucleons are coupled to finite-range me- 
son fields: isoscalar scalar, vector, isovector vector, and 
the Coulomb field and where the density dependence is 
modeled only by non-linear couplings of the scalar field 
[29l |30T | . This "standard" model is too constrained in the 
isovector channel and with respect to effective mass. It 
produces too narrow covariance ellipsoids for our correla- 
tion analysis. Therefore we augmented it by tensor cou- 
plings of vector fields [2!| and by an isovector scalar field 
with mass 980 MeV, denoting the resulting functional as 
RMF-5-t. We fit the model parameters to the same pool 
of data as SV-min. Since the resulting NMP of RMF-6- 
t (#=197MeV, a sym 38MeV, m*/m=0.59) strongly de- 
viate from the accepted values (as all traditional RMF 
models) we use this model only to discuss the robustness 
of our certain predictions and to illustrate the model de- 
pendence of the statistical analysis. 

Results. In our study, we studied selected NMP and 
a number of observables related to isovector properties 
of finite nuclei such as neutron skins and radii, binding 
energy differences, and dipole polarizability [31(. The 
latter one, the key quantity for static response, has been 
calculated within the RPA method: 



a D = 2 {\{®n\b\®a)\ 2 /E n ), (7) 

nGRPA 



where n runs over the excitation spectrum, E n is the 
excitation energy of the RPA state and D the elec- 

tric dipole operator (see Ref. [7[ for details of RPA 
calculations). We also investigated the energies of Gi- 
ant Resonances: monopole (GMR), dipole (GDR), and 
quadrupole (GQR), and the low-energy dipole strength 
in neutron-rich nuclei: 



B{E1; PDR) = Y B(El,ri) 
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with £" max =10MeV. The latter quantity is sometimes re- 
lated to the PDR strength pl420j. 

Figure [T] shows covariance ellipsoids for two pairs of 
observables in 208 Pb that nicely illustrate the cases of 
strong correlation (r s ki n and an; cab=0.98) and weak 
correlation (r s ki n and m*/m; cab=0.11). Figure [5] shows 
correlations with the point-neutron distribution form fac- 
tor F n in 208 Pb at (7=0.45fm _1 corresponding to PREX 
measurement. As expected, F n is strongly correlated 
with r s kinj r n mS ' as wen as with neutron skins in other 
neutron-rich nuclei. Almost equally strong is the cor- 
relation with the dipole polarizability. Not surprisingly, 
one can see excellent correlation of F n with NMP: a sym , 
Q,g ym , and d(E/A) n /dp n . All those quantities can thus 
be viewed as good indicators of isovector properties of 
nuclei. 

The excellent correlation between the neutron skin 
and dipole polarizability is not surprising as r s ki n oc 
otDdsym [HI- The experimental value of ao for 208 Pb is 
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FIG. 1: (Color online) The covariance ellipsoids for two pairs 
of observables as indicated. The filled area shows the region of 
reasonable domain p. Left: neutron skin and isovector dipole 
polarizability in 208 Pb. Right: neutron skin in 208 Pb and 
effective nucleon mass m*/m in symmetric nuclear matter. 
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FIG. 2: (Color online) Correlation ((6]) of various observables 
with the neutron form factor F n (q=0A5 fm _1 ) in 208 Pb. 



13.3±1.4fm 2 /MeV [H while the value obtained by the 
Lorentz fit to the total experimental photodisintcgration 
cross section is I3.6fm 2 /MeV [3lJ. As seen in Fig. [TJ 
both values are consistent with the SV-min predictions 
for r s kin- However, a 10% experimental uncertainty due 
to statistical and photon-beam calibration errors makes 
it impossible to use the current best value of old as an 
independent check on neutron skin. 
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The nuclear and neutron matter binding energy seem 
poorly correlated with F n , in accordance with Ref. 0. 
Our covariance analysis suggests the lack of correlation 
between F n (or neutron skin) and PDR strength; GMR, 
GDR, and GQR energies; isoscalar and isovector effective 
mass, incompressibility, and saturation density (see also 
Refs. [H El)- Those quantities can thus be viewed as 
poor isovector indicators. 

According to calculations, the degree of correlation 
with F n in 208 Pb in general deteriorates with decreasing 
mass number (see also discussion in [5l 111]). This can be 
explained in terms of increased importance of shell effects 
in lighter nuclei. Shell effects are also responsible for the 
lack of correlation between r s ki n and PDR strength. The 
low-energy El strength is greatly impacted by the de- 
tailed single-particle structure around the Fermi level and 
thus varies rapidly with EDF parameters. This reduces 
correlation with quantities which are weakly influenced 
by shell effects, e.g., NMP. The results of correlations 
for RMF-(5-t are very similar to those from SV-min. In 
particular, large correlations with neutron radii are pre- 
dicted for a sym , and a' sym and small for E/A, K, and 
m* jm. 

To estimate the impact of precise experimental deter- 
mination of neutron skin, we generated a new functional 
SV-min-i?„ by adding the value of neutron radius in 
208 Pb, r™ s =5.61 fm, with an adopted error 0.02 fm (0.4% 
measurement) and 0.05 fm (1% measurement), to the set 
of fit observables. (The main difference between SV-min 
and SV-min- R n is a slight reduction of isovector NMPs in 
the latter EDF: a sym from 30.66 MeV to 30.51 MeV; a' sym 
from 92.73 MeV fm 3 to 89.85 MeVfm 3 ; and k T rk from 
0.0765 to 0.057.) Assuming a 0.4% uncertainty in r™ ls , 
calculated uncertainties on isovector indicators shrink by 
about a factor of two. Figure [3] illustrates this tendency: 
it compares extrapolation errors for the neutron matter 
EOS in EDF SV-min and SV-min-i?„. The impact of a 
1% measurement is much smaller, at least for the range 
of densities considered. 

We also carried out calculations with a new EDF ob- 
tained by a new fit where the neutron-rich nuclei have 
been given more weight (a factor 2 to 3 for the three 
outermost neutron rich isotopes in most chains). The 
purpose of this exercise is to simulate the expected in- 
creased amount of data on neutron-rich nuclei. While the 
correlations seem to change very little, the extrapolation 
uncertainties in neutron observables shrink by a factor 
of 1.5-2.0. For instance, with this new functional, the 
predicted neutron skin in 208 Pb is r s ki n =0.f 91(0.024) fm, 
as compared to SV-min value of r s ki n =0. 170(0.037) fm. 
This exercise demonstrates that detailed conclusions of 
the statistical analysis depend on a chosen model and a 
selected set of fit observables. This point is also illus- 
trated in Fig. [3] the neutron matter EOS predicted in 
SV-min is very different from that of RMF-<5-t (which, as 
discussed earlier, yields unphysical NMP). 

Conclusions. In summary, we propose to use a sta- 
tistical least-squares analysis to identify the impact of 
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FIG. 3: (Color online) Extrapolation errors for the neutron 
matter EOS predicted by EDF SV-min (obtained by a fit 
to the standard pool of data) and SV-min-i?„ (obtained by 
adding to the data set the neutron radius in 208 Pb with an 
adopted error of 0.02 fm and 0.05 fm. The neutron EOS pre- 
dicted by RMF-<5-t is also shown for comparison). 



new observables, quantify correlations between predicted 
observables, and assess uncertainties of theoretical pre- 
dictions. To illustrate the concept, we studied the the 
neutron radius of of 208 Pb. By means of covariance anal- 
ysis we identified a set of good isovector indicators that 
correlate very well with the neutron form factor of 208 Pb. 
These are: neutron skins and radii in neutron-rich nuclei, 
dipole polarizability, and the nuclear matter properties 
such as symmetry energy and pressure. An indicator 
that is particularly attractive, as it can be measured in 
finite nuclei, is dipole polarizability. Unfortunately, the 
current best experimental value of qd in 208 Pb is not 
known precisely enough to offer an independent check on 
neutron skin or to provide a quality constraint on EDF. 
We also demonstrate that nuclear and neutron matter 
binding energy, low-energy El strength, giant resonance 
energies, isoscalar and isovector effective mass, incom- 
pressibility, and saturation density are poor indicators of 
isovector properties, at least those related to r s ki n - 

We discussed the impact of PREX measurement on 
theoretical uncertainties for neutron-rich nuclei or neu- 
tron matter and concluded that it will provide a valu- 
able constraint on the nuclear energy functional that will 
reduce theoretical error bars on the neutron-rich side. 
While we have good reason to believe that our general 
conclusion about good and poor isovector indicators is 
robust, predictions for individual observables are obvi- 
ously model dependent, as shown in Fig. |3] This is an 
important point: even the best statistical analysis is not 
going to eliminate systematic errors due to incorrect the- 
oretical assumptions. 

While our discussion is pertaining to the nuclear DFT, 
as the DFT is an obvious tool of choice to handle com- 
plex heavy nuclei and neutron skins, we believe that the 
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methodology used in this work should be of interest to 
any theoretical framework that contains parameters fine- 
tuned to experiment. Examples include fits of nucleon- 
nucleon forces to scattering and few-body data, adjust- 
ments of shell-model matrix elements, fits of coupling 
constants of symmetry-dictated Hamiltonians. 
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